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ABSTRACT 

We present here additional results of a spectropolarimetric survey of a small sample 
of stars ranging from spectral type MO to M8 aimed at investigating observationally 
how dynamo processes operate in stars on both sides of the full convection threshold 
(spectral type M4). 

The present paper focuses on early M stars (M0-M3), i.e. above the full convec- 
tion threshold. Applying tomographic imaging techniques to time series of rotationally 
modulated circularly polarised profiles collected with the NARVAL spectropolarime- 
ter, we determine the rotation period and reconstruct the large-scale magnetic topolo- 
gies of 6 early M dwarfs. We find that early-M stars preferentially host large-scale 
fields with dominantly toroidal and non-axisymmetric poloidal configurations, along 
with significant differential rotation (and long-term variability); only the lowest-mass 
star of our subsample is found to host an almost fully poloidal, mainly axisymmetric 
large-scale field ressembling those found in mid-M dwarfs. 

This abrupt change in the large-scale magnetic topologies of M dwarfs (occuring at 
spectral type M3) has no related signature on X-ray luminosities (measuring the total 
amount of magnetic fiux); it thus suggests that underlying dynamo processes become 
more efficient at producing large-scale fields (despite producing the same fiux) at 
spectral types later than M3. We suspect that this change relates to the rapid decrease 
in the radiative cores of low-mass stars and to the simultaneous sharp increase of the 
convective turnover times (with decreasing stellar mass) that models predict to occur 
at M3; it may also be (at least partly) responsible for the reduced magnetic braking 
reported for fully-convective stars. 

Key words: stars: magnetic fields - stars: low-mass - stars: rotation - stars: activity 
- techniques: spectropolarimetry 
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1 INTRODUCTION 

Activity and magnetic fi elds are ubiquitous t o cool stars of 
all sp ectral types (e.g., ISaar Linskvl [TqSSI : iDonati et al.l 
119971 ). The analogy with the Sun suggests that these fields 
are produced through dynamo mechanisms operating in a 
thin interface layer at the base of the convective envelope, 
where angular rotation gradients are strongest. Starting 
from a weak poloidal configuration, differential rotation pro- 
gressively winds the field around the star, building a strong 
toroidal belt at the base of the convective zone; cyclonic tur- 
bulence then restores a weak poloidal field (with a polarity 
opposite to that of the initial o ne) once the toroidal field has 
grown unstable (|Parker|[T955h . However, there is still con- 
siderable controversy (even for the Sun itself) on how and 
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where exactly the field is amplified and on which physical 
processes (differential rotation, meridional ci rculation) are 
mainly controlling the magnetic cycle (e.g., ICharbonneaul 

[iooi). 

Observing stars other than the Sun is of obvious inter- 
est for this question as they provide a direct way of study- 
ing how dynamo processes depend on fundamental stellar 
parameters such as mass (and thus convective depth) and 
rotation rate. In this respect, low-mass fully- convective stars 
are particularly interesting as they host no interface layer 
(where dynamo processes presumably concentrate in the 
Sun), and nevertheless show bo th strong magnetic fields 
fe.g.Jjohns-Krull Vale^ ll996) and intense activity (e.g., 
iDelfosse et al.l 19981 ). Manv theoretical mod els were pro- 
posed to attempt resolving this issue (e .g., iDurnev et al.l 
I1993I : iDobler et al.| [2006l : iBrownind I2OO8I ) but observations 
of the large-scale magnetic fields at the surfaces of fully- 
convective stars, and in particular of their poloidal and 
toroidal components, are still rare. 

Thanks to time-resolved spectropolarimetric observa- 
tions of cool stars, we are now able to recover informa- 
tion on how magnetic fields distribute at (and emer ge from) 
the surfaces of cool active stars fe.g., IPonati et al.ll2003al ). 
With the advent of new generation instruments optimised 
in this very purpose (ESPaDOnS at the 3.6m Canada- 
France-Hawaii Telescope and N ARVAL at the 2m Telescope 
Bernard Lyot in France, iDona ti 2003), this technique can 
now be also applied to faint M dwarfs. The first attempt 
focussed on the fully- convective rapidly-rotating M4 star 
V374 Peg and revealed that, against all theoretical expec- 
tations, fully-convective M dwarfs are apparently very effi- 
cient at produ cing strong large-sc ale mainly-axisymmetric 
poloidal fields ([Ponati et al.ll2006a ) despite very small lev- 
els of differential rotation. A follow-up study confirmed this 
point and further demonstrated that the magnetic config- 
uration of V374 Peg is ap parently stable on timescales of 
:^ 1 yr flMorin et al.il2008al ). 

To investigate this issue in more details, we embarked in 
a spectropolarimetric survey of a small sample of M dwarfs, 
located both above and below the full-convection thresh- 
old (corresponding to a mass of 0.35 M© and to a spec- 
tral type of M4). The first results of this survey, focussing 
mainly on active M4 dwarfs, confirms the results obtained on 
V374 Peg and demonstrate that strong large-scale mainly- 
axisymmetric poloidal fields are indeed fairl y common in 
fully-convective mid-M active dwarfs (iMorin et al.l l2008bl , 
hereafter M08). In this new paper, we present the results for 
the 6 early M dwarfs that we observed in this survey, namely 
DT Vir, DS Leo, CE Boo, OT Ser, GJ 182 and GJ 49, with 
spectral types ranging from MO to M3 (see Table [1]). After 
briefly describing the observations and the magnetic mod- 
elling method that we use, we detail the results obtained for 
each star and discuss their implication for our understanding 
of dynamo processes in cool active stars. 



2 OBSERVATIONS 

Spectropolarimetric observations of the selected M dwarfs 
were collected with NARVAL and the 2m Telescope Bernard 
Lyot (TBL), between 2007 Jan and 2008 Feb (in 3 dif- 
ferent runs). NARVAL is a twin spectropolarimeter copied 
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Figure 1. Average unpolarised LSD profiles for DT Vir (red line), 
DS Leo (green), CE Boo (dark blue), OT Ser (light blue), GJ 182 
(pink) and GJ 49 (yellow). All profiles are centred on the rest 
velocity. 



from ESPaDOnS ('Donati"2003') , yielding full coverage of 
the optical domain (370 to 1000 nm) at a resolving power of 
about 65 000 in a single exposure. Each polarisation expo- 
sure consists of 4 individual subexposures taken in different 
polarimeter configurations and combined together to filter 
out all spurious p olarisation signatures at first order (e.g., 
iDonati et al]ll997l ). 

Data reduction was carried out using Libre ESpRIT, 
a fully automated package/pipeline installed at TBL and 
performing optimal extraction of unpolarised (Stokes /) 
and circularly polar ised (Stokes V) spectra as described in 
IDonati et al T (|l997l ). The peak signal-to-noise ratios (S/N) 
per 2.6 kms~^ velocity bin that we obtained in the collected 
spectra range from 100 to 400 depending on the magnitude 
and weather conditions. The journal of observations for all 
stars is presented in Tables [2] to [71 

All spectra are automatically corrected from spectral 
shifts resulting from instrumental effects (eg mechanical flex- 
ures, temperature or pressure variations) using telluric lines 
as a reference. Though not perfect, this procedure allows 
spectra to be secured with a radi al velocity (RV) pre cision 
of better than 0.030 kms"^ fe.g.. iMoutou et al.l[2007l) 

Least-squares deconvolution (LSD, iDonati et al.lll997l ) 
was applied to all spectra to extract the polarisation signal 
from most photospheric atomic lines and compute a mean 
Zeeman signature corresponding to an average photospheric 
profile (centred at 700 nm and with an effective Lande fac- 
tor of 1.2). The line list used in this process is derived from 
an A tlas9 local thermodynamic equilibrium model (|Kuruc3 
I1993I ) matching the properties of our sample, and includes 
about 5,000 moderate to strong atomic lines (i.e., with a rel- 
ative depth larger than 40% prior to any macroscopic broad- 
ening). The resulting multiplex gain in S/N is about 15 (see 
Table [2]). Zeeman signatures are detected in most cases on 
all stars. Radial velocities are obtained from Gaussian fits 
to all LSD unpolarised profiles. Average unpolarised LSD 
profiles of all stars are shown in Fig. [1] 
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Table 1. Fundamental parameters of our sample early-M dwarf stars. Columns 1 to 6 respectively list the name, the spectral type, the 
mass (derived from the Hipparcos distance and the J, H and K magnitudes using the mass-luminosity relations of Delfosse e t al.l 
[2000 except for GJ 182; for which we used the evolutiona ry models of Baraffe et al. 1998, see text), the logarithmic bolometric lumi- 
nosity logLboi (de rived from the mass a nd the models of iBaraffe et al . 1998), the logarithmic relative X-ray luminosity log i?x (i-e., 
log Lx/Lbol, from lKiraga Stepienll2007l or from the NEXXUS data base, iSchmitt Liefke 2004, or from Woo d et al]|l994l for GJ 49), 
the projected rotation v elocity (this paper, accuracy ~ 1 kms~^), the rotation period (this paper), the convective turnover time Tq (from 
Kiraga SteDienll2007h . the effective Rossby number Ro = Prot/Tc, the radius Ri, (predicted by the theoretical models of iBaraffe et al.l 
19981 ) and the assumed inclination angle of the rotation axis to the line-of-sight (this paper). 



Star ST log Lboi log v sin i Prot Tc Ro R-^ i 

(Me) (ergs-i) (kms-i) (d) (d) (Rq) (°) 



GJ 182 


M0.5 


0.75 


32.7 


-3.1 


10 


4.35 


25 


0.174 


0.82 


60 


DT Vir / GJ 494A 


MO. 5 


0.59 


32.3 


-3.4 


11 


2.85 


31 


0.092 


0.53 


60 


DS Leo / GJ 410 


MO 


0.58 


32.3 


-4.0 


2 


14.0 


32 


0.438 


0.52 


60 


GJ 49 


Ml. 5 


0.57 


32.3 


< -4.3 


1 


18.6 


33 


0.564 


0.51 


45 


OT Ser / GJ 9520 


Ml. 5 


0.55 


32.2 


-3.4 


6 


3.40 


35 


0.097 


0.49 


45 


CE Boo / GJ 569A 


M2.5 


0.48 


32.1 


-3.7 


1 


14.7 


42 


0.350 


0.43 


45 



Table 2. Journal of observations for DT Vir=GJ 494A. Columns 1-7 list the UT date, the heliocentric Julian date, the UT time, 
the total exposure time, the peak signal to noise ratio (per 2.6 kms"-*^ velocity bin) and the rms noise level in the LSD Stokes V 
profile (relative to the unpolarised continuum level and per 1.8 kms"-*^ velocity bin). In col. 8 and 9, we list the rotational cycles (using 
ephemeris HJD= 2,454,100.0 + 2.85E) as well as the longitudinal fields (with error bars) and the radial velocities Vr (absolute 
accuracy 0.10 kms"-*^, internal accuracy 0.03 kms"-*^) associated to each exposure. 



Date 


HJD 


UT 


^exp 


S/N 


C^LSD 


Cycle 


Be 








(2,454,000+) 


(h:m:s) 


(s) 




(10-4/e) 




(G) 


(kms-1) 


2007 Jan 26 


126.70173 


04:47:56 


4x800 


240 


3.3 


9.369 


-17.6 ±9.3 


-13.32 


Jan 27 


127.67710 


04:12:22 


4x800 


310 


2.3 


9.711 


-3.8 ±6.7 


-13.44 


Jan 28 


128.69001 


04:30:50 


4x800 


350 


2.0 


10.067 


66.0 ±5.8 


-13.04 


Jan 29 


129.65374 


03:38:29 


4x800 


340 


2.1 


10.405 


-27.6 ±6.0 


-13.29 


Feb 02 


133.71569 


05:07:13 


4x800 


350 


2.0 


11.830 


23.3 ±5.7 


-13.29 


Feb 03 


134.69961 


04:43:57 


4x800 


340 


2.1 


12.175 


49.7 ±6.0 


-13.00 


Feb 04 


135.71225 


05:02:02 


4x800 


360 


1.9 


12.531 


-31.1 ±5.6 


-13.32 


Feb 05 


136.66208 


03:49:41 


2x800 


140 


5.2 


12.864 


28.0 ±15.5 


-13.25 


2007 Dec 


28 


462.71536 


05:11:18 


4x750 


330 


2.2 


127.269 


30.5 ±6.5 


-13.64 


Dec 


29 


463.74605 


05:55:22 


4x750 


360 


2.0 


127.630 


15.5 ±5.7 


-13.67 


Dec 


31 


465.74401 


05:52:10 


4x700 


290 


2.5 


128.331 


29.0 ±7.1 


-13.66 


2008 Jan 


01 


466.74583 


05:54:39 


4x700 


360 


2.0 


128.683 


24.5 ±5.7 


-13.69 


Jan 


19 


484.63507 


03:12:52 


4x600 


170 


4.5 


134.960 


26.3 ± 12.8 


-13.57 


Jan 


20 


485.61280 


02:40:40 


4x600 


200 


4.0 


135.303 


29.4 ± 11.3 


-13.77 


Jan 


23 


488.62610 


02:59:27 


4x600 


250 


3.1 


136.360 


23.6 ±9.0 


-13.74 


Jan 


24 


489.62041 


02:51:08 


4x600 


240 


2.9 


136.709 


1.5 ±8.6 


-13.65 


Jan 


26 


491.60391 


02:27:08 


4x600 


230 


3.3 


137.405 


39.8 ±9.5 


-13.69 


Jan 


27 


492.62606 


02:58:54 


4x600 


320 


2.2 


137.763 


-4.5 ±6.2 


-13.67 


Jan 


28 


493.63758 


03:15:22 


4x600 


300 


2.4 


138.118 


4.5 ±6.7 


-13.56 


Feb 


03 


499.63751 


03:14:34 


4x600 


280 


2.7 


140.224 


8.9 ±7.7 


-13.80 


Feb 


05 


501.64773 


03:29:03 


4x600 


290 


2.4 


140.929 


53.9 ±7.0 


-13.59 


Feb 


06 


502.63636 


03:12:34 


4x600 


320 


2.2 


141.276 


-7.1 ±6.4 


-13.75 


Feb 


07 


503.64038 


03:18:15 


4x600 


290 


2.4 


141.628 


23.2 ±7.1 


-13.69 


Feb 


10 


506.64784 


03:28:40 


4x600 


310 


2.3 


142.684 


7.0 ±6.6 


-13.65 


Feb 


12 


508.63939 


03:16:18 


4x600 


300 


2.4 


143.382 


-4.9 ±6.9 


-13.82 


Feb 


14 


510.63942 


03:16:08 


4x600 


230 


3.2 


144.084 


7.4 ±9.2 


-13.61 


Feb 


15 


511.65542 


03:39:04 


4x600 


280 


2.5 


144.440 


6.7 ±7.3 


-13.77 


Feb 


16 


512.64053 


03:17:32 


4x600 


300 


2.4 


144.786 


-2.0 ±6.7 


-13.68 


Feb 


17 


513.64033 


03:17:09 


4x600 


310 


2.3 


145.137 


-0.1 ±6.6 


-13.65 



3 MAGNETIC MODELLING 

The magnetic model we use to describe the time series of 
rotationally modulated LSD Stokes V profiles is described 
in MOB. We recall it briefly here and refer the reader to MOB 
for further details. 



To des cribe the magnetic fi eld, we use the same de- 
scription as iDonati et al.l (|2006bh . The field is decomposed 
into its poloidal and toroidal components, both expressed 
as spherical harmonics expansions. The imaging process is 
based on the principles of maximum entropy image recon- 
struction, with entropy (i.e., quantifying the amount of re- 
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Table 3. Same as Table[2]for DS Leo=GJ 410. Rotational cycles are now computed using ephemeris HJD= 2,454, 100.0 + 14. 0£^. 



Date 


HJD 


UT 


t'exp 


S/N 




Cycle 












(h.m.s) 






J-c ) 






(km s"-*-) 


2007 Jan 26 


126.65015 


03:31:14 


4x800 


250 


2.7 


1.904 


-33.9 ±4.1 


-14.08 


Jan 27 


127.63491 


03:09:13 


4x800 


270 


2.3 


1.974 


-15.9 ±3.4 


-14.05 


Jan 28 


128.64692 


03:26:27 


4x800 


320 


2.0 


2.046 


-1.3 ±3.0 


-14.08 


Jan 29 


129.61070 


02:34:13 


4x800 


330 


1.9 


2.115 


13.1 ±2.9 


-14.06 


Jan 30 


130.65362 


03:35:57 


4x800 


340 


1.7 


2.189 


29.7 ±2.6 


-14.08 


Feb 02 


133.67041 


03:59:57 


4x800 


340 


1.8 


2.405 


6.3 ±2.8 


-14.03 


Feb 03 


134.65677 


03:40:15 


4x800 


340 


1.8 


2.475 


7.0 ±2.8 


-14.01 


Feb 04 


135.66947 


03:58:28 


4x800 


350 


1.8 


2.548 


3.7 ±2.7 


-14.02 


Feb 05 


136.63031 


03:02:02 


4x800 


200 


3.2 


2.617 


2.4 ±5.1 


-14.04 


2007 Dec 


28 


462.67598 


04:11:04 


4x750 


320 


2.0 


25.905 


9.9 ± 2.9 


-14.09 


Dec 


29 


463.70487 


04:52:33 


4x750 


340 


1.8 


25.979 


1.5 ± 2.7 


-14.09 


Dec 


31 


465.70955 


04:59:06 


4x600 


260 


2.5 


26.122 


8.6 ± 3.7 


-14.08 


2008 Jan 


01 


466.70912 


04:58:22 


4x700 


350 


1.8 


26.194 


24.5 ± 2.7 


-14.09 


Jan 


02 


467.68306 


04:20:44 


2x600 


130 


5.1 


26.263 


13.9 ± 7.7 


-14.07 


Jan 


03 


468.70748 


04:55:48 


4x700 


180 


3.8 


26.336 


26.9 ±5.6 


-14.05 


Jan 


07 


472.62578 


02:57:46 


4x700 


160 


4.3 


26.616 


-22.1 ± 7.0 


-14.06 


Jan 


08 


473.62319 


02:53:56 


4x700 


260 


2.4 


26.687 


-4.0 ± 3.9 


-14.11 


Jan 


19 


484.58050 


01:51:30 


4x600 


180 


4.0 


27.470 


-17.2 ± 5.9 


-14.04 


Jan 


20 


485.55637 


01:16:40 


4x600 


190 


3.7 


27.540 


-26.8 ± 5.5 


-14.05 


Jan 


22 


487.58395 


01:56:14 


4x600 


230 


2.8 


27.685 


-9.6 ±4.4 


-14.04 


Jan 


24 


489.53279 


00:42:24 


4x600 


250 


2.6 


27.824 


6.4 ±3.8 


-14.05 


Jan 


26 


491.54899 


01:05:35 


4x600 


240 


2.6 


27.968 


7.1 ± 3.8 


-14.08 


Jan 


27 


492.57173 


01:38:15 


4x600 


310 


2.0 


28.041 


-3.5 ± 3.0 


-14.07 


Jan 


28 


493.58339 


01:54:59 


4x600 


290 


2.2 


28.113 


6.7 ±3.3 


-14.04 


Feb 


03 


499.60160 


02:20:49 


4x600 


270 


2.3 


28.543 


-22.3 ±3.8 


-14.11 


Feb 


05 


501.58341 


01:54:31 


4x600 


250 


2.6 


28.684 


-0.9 ±4.1 


-14.08 


Feb 


06 


502.58201 


01:52:26 


4x600 


310 


2.1 


28.756 


4.3 ±3.1 


-14.08 


Feb 


07 


503.58042 


01:50:06 


4x600 


290 


2.2 


28.827 


3.5 ±3.2 


-14.10 


Feb 


10 


506.59331 


02:08:31 


4x600 


290 


2.1 


29.042 


4.2 ±3.2 


-14.03 


Feb 


11 


507.55697 


01:16:08 


4x600 


300 


2.0 


29.111 


5.8 ±2.9 


-14.06 


Feb 


12 


508.58557 


01:57:17 


4x600 


310 


2.2 


29.185 


16.1 ±3.3 


-14.04 


Feb 


13 


509.59071 


02:04:38 


4x600 


290 


2.2 


29.256 


16.1 ±3.3 


-14.02 


Feb 


14 


510.58572 


01:57:25 


4x600 


230 


2.9 


29.328 


10.9 ±4.3 


-14.03 


Feb 


15 


511.60258 


02:21:39 


4x600 


300 


2.1 


29.400 


-9.4 ±3.2 


-13.98 


Feb 


16 


512.58834 


02:01:06 


4x600 


270 


2.4 


29.471 


-16.5 ±3.5 


-14.04 


Feb 


17 


513.58759 


01:59:59 


4x600 


300 


2.1 


29.542 


-16.9 ±3.1 


-14.02 



constructed information) being calculated from the coeffi- 
cients of the spherical harmonics expansions. Starting from 
a null magnetic field, we iteratively improve our magnetic 
model by comparing the synthetic Stokes V profiles with 
the observed ones, until we reach an optimal field topology 
that reproduces the data at a given ^ level (i.e., usually 
down to noise level, corresponding to a unit reduced level 
= 1). The inversion problem being partly ill-posed, we 
use the entropy function to select the magnetic field with 
lowest information content among all those reproducing the 
data equally well. Given that most stars considered here ro- 
tate no more than moderately, we limit spherical harmonics 
expansions to ^ < 10, usually up to 8 for moderate rotators 
(e.g., DT Vir) and up to 5 only for the slower ones. In all 
cases (even the slowest rotators), we need to set > 3 to 
reproduce successfully the data at noise level. 

To compute the synthetic profiles corresponding to a 
given magnetic topology, we divide the surface of the star 
into a grid of elementary surface cells (typically 5,000), in 
which the 3 components of the magnetic field (in spherical 
coordinates) are estimated directly from the spherical har- 
monics expansions used to describe the field. Using Unno- 



Rachkovsky's equations (e.g., iLandi degl'Innocentil [iQQ^ l , 
we compute the contribution of each grid cell to the Stokes 
V profiles and integrate all contributions from the visible 
stellar hemisphere at each observed rotation phase. The free 
parameters in Unno-Rachkovsky's equations (describing the 
shape of the unpolarised line profile from a non-magnetic 
grid cell) are obtained by fitting the Stokes / LSD profiles 
of a very- slowly rotating and weakly active star of similar 
spectral type (e.g., GJ 205). 

Reproducing both the amplitude and shape of LSD 
Stokes V profiles in the particular case of stars with strong 
fields and sharp lines requires that we introduce a filling 
factor (called /) describing the fractional amount of flux 
producing circular polarisation (assumed constant over the 
whole star). At first glance, this may seem in contradic- 
tion with the fact that we are mostly sensitive to large-scale 
fields, i.e., fields whose spatial coherency is much larger than 
the size of our grid cells; however, large-scale fields can po- 
tentially be also structured on a small-scale, e.g., with con- 
vection compressing the field into a small section of each 
cell but keeping the flux constant over the cell surface. We 
suspect that this is the case here. In practice, it means that 
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Table 4. Same as Table [2] for CE Boo=GJ 569A. Rotational cycles are now computed using ephemeris HJD= 2,454, 100.0 + \A.7E. 



Date 


HJD 


UT 


^exp 


S/N 


C^LSD 


Cycle 


Be 


Vr 




(2,454,000+) 


(h:m:s) 


(s) 




(10-4/e) 




(G) 


(kms-1) 


9nnR Tan 1Q 


^O^.UUU^i- 


flzl-ni -AA 


4x600 


130 


6.2 


26.168 


_Qq _|_ 1 4 c: 


—7.31 


Tan 9n 


AR^ (\RO^R 


r\A-9A-^A 


4x600 


200 


3.9 


26.237 




—7.36 


Tan 9*^ 


ARR fit=;s9n 


n'^-4Q-97 

UO.^c/. Z ( 


4x600 


230 


3.3 


26.439 


_c:q 0170 


—7.28 


Tan 94 
J dii 


ARQ fit=;i Q7 


UO.^U.Zi. 


4x600 


190 


3.7 


26.507 


_7q q _|_ 
— ( 0.0 ziz 0.0 


—7.29 


Tan 9fi 
«j dii ZiU 


4Q1 fiSfifi7 




4x600 


170 


4.7 


26.642 


_i 1 n _L 1 1 

_l_ J. J. . 


—7.31 


Tan 97 
J dii Z ( 


4Q9 fi^Sn9 


r\'^-AR-A9 


4x600 


270 


2.6 


26.711 


— c/\j.o zn u.i- 


—7.34 


Jan 28 


493.67096 


04:07:12 


4x600 


280 


2.5 


26.780 


-114.2 ±6.1 


-7.31 


Jan 30 


495.68078 


04:21:06 


4x600 


250 


2.9 


26.917 


-111.0 ±6.8 


-7.35 


Feb 03 


499.69324 


04:38:34 


4x600 


260 


2.7 


27.190 


-70.4 ±6.5 


-7.42 


Feb 05 


501.68094 


04:20:37 


4x600 


250 


2.9 


27.325 


-53.8 ±6.9 


-7.35 


Feb 07 


503.67364 


04:09:51 


4x600 


260 


2.7 


27.461 


-54.9 ±6.4 


-7.35 


Feb 10 


506.67974 


04:18:16 


4x600 


270 


2.7 


27.665 


-91.0 ±6.3 


-7.32 


Feb 11 


507.68809 


04:30:10 


4x600 


280 


2.4 


27.734 


-113.4 ±5.9 


-7.36 


Feb 12 


508.67116 


04:05:40 


4x600 


240 


3.0 


27.801 


-115.0 ±7.1 


-7.34 


Feb 13 


509.67882 


04:16:35 


4x600 


230 


3.1 


27.869 


-108.6 ±7.4 


-7.33 


Feb 14 


510.67193 


04:06:32 


4x600 


200 


3.7 


27.937 


-106.7 ±8.6 


-7.35 


Feb 15 


511.68789 


04:29:24 


4x600 


270 


2.7 


28.006 


-111.7±6.4 


-7.31 


Feb 16 


512.67287 


04:07:39 


4x600 


260 


2.6 


28.073 


-85.5 ±6.2 


-7.37 


Feb 17 


513.67252 


04:07:02 


4x600 


270 


2.6 


28.141 


-79.9 ±6.2 


-7.34 



the circularly polarised flux that we get from each grid cell 
is given by fV\oc where Vioc is the Stokes V profile derived 
from Unno-Rachkovsky's equations for a magnetic strength 
of B/ f. In practice, this simple model ensures that both the 
width and amplitude of Stokes V signatures can be fitted 
simultaneously; typical values of / range from 0.10 to 0.15 
in active mid-M dwarfs with sharp lines (MOB). For stars ro- 
tating more rapidly and/or hosting intrinsically weak fields, 
different values of / produce very similar fits and undistin- 
guishable magnetic flux maps, in which case we arbitrarily 
set f = 1- This model has proved rather successful at repro- 
ducing the observed times series of Stokes V profiles in mid- 
M dw arfs (MOB) and classical T Tauri stars (Donati et al.l 
I2OOBI ): we therefore use it again for the present study. 

We can also implement differential rotation for com- 
puting the synthetic Stokes V profiles corresponding to our 
magnetic model. For this purpose, we use a Sun-like surface 
rotation pattern with the rotation rate varying with latitude 
as Qeq — dQsin^ 0, Qeq being the angular rotation rate at 
the equator and dQ the difference in angular rotation rate 
between the equator and the pole. By carrying out recon- 
structions (at constant information content) for a range of 
Qeq and dQ values, we can investigate how the fit quality 
varies with differential rotation; differential rotation is de- 
tected when the of the fit to the data shows a well defined 
minimum in the explored Qeq-o?^ domain, with the position 
of the minimum and the curvature of the surface at this 
point yield ing the optimal Qeg a nd dQ values and respective 
error bars (|Donati et al.ll2003bl l. 



4 DT VIR = GJ 494A = HIP 63942 

DT Vir is a magnetically active MO. 5 dwarf showing sig- 
nificant ro tational broadenin g in spectral lines {v sin i :^ 
10 kms~^, [Beuzit et al.ll2004l l: photometric variability sug- 
gests a short rotation period Prot of about 2.9 d in good 
agreement with its membership to the young galactic disc 



(iKiraga SteDienll2007h . The Hipparcos distance is 11.43 ± 
0.20 pc. It belongs to a distant binary system with an as- 
trometric period of about 14.5 yr (Heintz 1994) and a semi- 
major axis of 5-6 AU; the companion is about 4.4 magni- 
tudes fainter in K a nd is either a very- low- mass star or a 
young brown dwarf (|Beuzit et all |2004| . depending on the 
exact age of the system). U sing the mass- luminosity rela- 
tions of Delfosse et al. we estimate that the mass of 
DT Vir is 0.59 ± 0.02 M©. Given vsini and Prot, we infer 
that i^^sini is about 0.6 R©, i.e., already 10% larger than 
the radius expected from theoretical models (see Table [T]); 
we therefore set i = 60° in the imaging process (the result 
being weakly sensitive to variations of i of ±10°). 

Stokes V data were collected at 2 epochs, providing only 
moderate coverage of the rotation cycle at the first epoch 
but a dense and redundant coverage at the second epoch 
(see Table [2]). Stokes V signatures are clearly detected at all 
time, eventhough the corresponding longitudinal fields are 
usually low (ranging from -30 G to 70 G), much lower than 
those reported on mid-M dwarfs in particular (MOB). The 
projected rotation velocity that we derive from the Stokes 
/ profiles is v sin z = 11 ± 1 km s~^, in good agreement with 
the estimate of iBeuzit et al The RV we measure 

are different at both epochs, equal to —13.25 kms~^ and 
— 13.6B kms~^ respectively (with an absolute accuracy of 
about 0.10 kms~^); this difference likely refiects the binary 
motion. Moreover, the relative dispersion about the mean 
RV, respectively equal to 0.14 kms~^ and 0.07 kms~^, is 
larger than the internal RV accuracy of NARVAL (about 
0.03 kms"\ e.g., Moutou et al.ll20Q7l ) and likely refiects the 
intrinsic activity RV jitter of DT Vir; in 2007 (i.e., when the 
internal RV dispersion is largest), we find that RVs corre- 
late reasonably well with longitudinal fields, suggesting that 
the RV fiuctuations are indeed due to the magnetic activity. 
Given the moderate strength of the field and the significant 
rotational broadening of DT Vir, there is no need of adjust- 
ing / to optimise the fit quality. 

Assuming solid body rotation, we find that the rotation 
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Table 5. Same as Table[2]for OT Ser=GJ 9520. Rotational cycles are now computed using ephemeris HJD= 2,454, 100.0 + 3.40^;. 



Date 


HJD 


UT 


^exp 


S/N 




Cycle 










(2,454,000+) 


(h:m:s) 


(s) 




(10-4/c) 




(G) 


(kms~^) 


2007 Jul 


26 


308.42253 


22:05:53 


4x900 


270 


2.7 


61.301 


80.0 ±5.7 


6.89 


Jul 


28 


310.38549 


21:12:45 


4x900 


340 


2.0 


61.878 


55.0 ±4.6 


6.75 


Jul 


29 


311.43298 


22:21:15 


4x900 


370 


1.8 


62.186 


73.4 ±4.1 


6.87 


Jul 


30 


312.41904 


22:01:17 


4x900 


380 


1.8 


62.476 


74.9 ±3.8 


6.82 


Jul 


31 


313.41995 


22:02:42 


4x900 


230 


3.2 


62.771 


67.0 ±7.0 


6.82 


Aug 


02 


315.42156 


22:05:14 


4x900 


350 


1.9 


63.359 


80.9 ±4.2 


6.90 


Aug 


03 


316.42155 


22:05:20 


4x900 


350 


1.9 


63.653 


67.8 ±4.2 


6.84 


Aug 


04 


317.42236 


22:06:36 


4x900 


230 


3.1 


63.948 


55.5 ±7.1 


6.77 


Aug 


09 


322.41667 


21:58:58 


4x900 


330 


2.1 


65.417 


77.7 ±4.5 


6.87 


Aug 


10 


323.41738 


22:00:06 


4x900 


250 


2.8 


65.711 


70.5 ±6.1 


6.82 


Aug 


14 


327.35676 


20:33:15 


4x900 


270 


2.7 


66.870 


54.5 ±6.0 


6.79 


2008 Jan 


19 


484.74579 


05:56:56 


4x600 


150 


5.2 


113.160 


52.7 ± 10.6 


6.78 


Jan 


20 


485.71587 


05:13:45 


4x600 


190 


3.9 


113.446 


1.4 ±8.1 


6.80 


Jan 


23 


488.69155 


04:38:25 


4x600 


180 


4.1 


114.321 


8.9 ±8.7 


6.81 


Jan 


24 


489.68897 


04:34:36 


4x600 


160 


4.5 


114.614 


69.1 ±9.6 


6.76 


Jan 


26 


491.67101 


04:08:31 


4x600 


200 


3.6 


115.197 


47.5 ±7.6 


6.78 


Jan 


27 


492.69173 


04:38:14 


4x600 


260 


2.8 


115.498 


22.4 ±5.8 


6.79 


Jan 


28 


493.70473 


04:56:51 


4x600 


280 


2.4 


115.796 


90.9 ±5.3 


6.81 


Jan 


30 


495.71492 


05:11:18 


4x600 


250 


2.8 


116.387 


-6.6 ±5.9 


6.80 


Feb 


03 


499.72743 


05:28:51 


4x600 


240 


2.9 


117.567 


47.7 ±6.2 


6.72 


Feb 


05 


501.71560 


05:11:36 


4x600 


260 


2.7 


118.152 


86.3 ±5.7 


6.71 


Feb 


07 


503.70655 


04:58:20 


4x600 


280 


2.5 


118.737 


79.7 ±5.4 


6.75 


Feb 


10 


506.71299 


05:07:16 


4x600 


260 


2.7 


119.621 


52.6 ±5.6 


6.77 


Feb 


11 


507.72197 


05:20:05 


4x600 


280 


2.5 


119.918 


108.3 ±5.4 


6.71 


Feb 


12 


508.70420 


04:54:23 


4x600 


240 


3.0 


120.207 


56.9 ±6.2 


6.75 


Feb 


13 


509.62257 


02:56:43 


4x600 


230 


3.1 


120.477 


9.8 ±6.6 


6.78 


Feb 


13 


509.75390 


06:05:49 


4x600 


270 


2.6 


120.516 


27.8 ±5.5 


6.77 


Feb 


14 


510.70531 


04:55:45 


4x600 


220 


3.3 


120.796 


92.5 ±6.8 


6.78 


Feb 


15 


511.72139 


05:18:47 


4x600 


270 


2.5 


121.094 


104.7 ±5.4 


6.73 


Feb 


16 


512.70638 


04:57:03 


4x600 


280 


2.5 


121.384 


10.3 ±5.4 


6.74 


Feb 


17 


513.70577 


04:56:03 


4x600 


270 


2.5 


121.678 


76.0 ±5.4 


6.77 



Table 6. Same as Table[2]for GJ 182. Rotational cycles are now computed using ephemeris HJD= 2,454, 100.0 ± 4.35£^. 



Date 


HJD 


UT 


^exp 


S/N 


C^LSD 


Cycle 








(2,454,000±) 


(h:m:s) 


(s) 




(10-4/e) 




(G) 


(kms-1) 


2007 Jan 21 


122.37018 


20:48:02 


4x900 


210 


3.2 


5.143 


-84.4 ±8.0 


19.23 


Jan 26 


127.28464 


18:45:22 


4x900 


190 


3.6 


6.272 


-40.6 ±9.2 


19.22 


Jan 27 


128.28237 


18:42:12 


4x900 


240 


2.6 


6.502 


-37.9 ±6.4 


19.59 


Jan 30 


131.27339 


18:29:36 


2x900 


55 


12.7 


7.189 


-65.3 ±36.4 


19.22 


Feb 01 


133.29669 


19:03:23 


4x900 


250 


2.4 


7.654 


18.6 ±5.9 


19.66 


Feb 03 


135.28638 


18:48:47 


4x900 


190 


3.4 


8.112 


-84.2 ±8.7 


19.21 


Feb 08 


140.29936 


19:08:04 


4x900 


250 


2.5 


9.264 


-42.6 ±6.5 


19.24 



period providing the best fit to the data is close to that de- 
rived from photometric variations but slightly different for 
each of the 2 data sets, about 2.90 d for the 2008 data and 
2.80 d for the 2007 data; we therefore selected Prot = 2.85 d 
as the mean rotation period with which we phased all data. 
Using this value of Prot yields a slightly chaotic phase de- 
pendence for the 2008 Bi data, with points at nearby phases 
but different cycles (e.g., on Dec 28 and Feb 06, i.e., at ro- 
tation cycles 127.269 and 141.276) showing discrepant field 
values; we suspect it indicates significant surface differen- 
tial rotation on DT Vir. Further confirmation comes from 
our finding that the 2008 data cannot be fitted down to 
— 1 for solid body rotation; using differential rotation, 



we are able to fit to the data down to noise level, with the 
surface (at given information content) showing a clear min- 
imum. The differential rotation parameters we obtain are 
Qeq = 2.200±0.003 radd-^ and dVt = 0.060±0.006 radd-\ 
corresponding to rotation periods at the equator and pole 
of 2.85 d and 2.94 d res pectively (bracketing the estimate of 
Kiraga & Stepienl l2007h . The photospheric shear of DT Vir 
is very similar to that of the Sun, with the equator lapping 
the pole by one rotation cycle every 105 ± 10 d. 

The optimal maximum entropy fit to the Stokes V data 
that we obtain (including the effect of differential rotation) 
is shown in Fig. [2] and corresponds to = 1? i-e., to a x? 
improvement over a non-magnetic model of x31 and x22 for 
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Figure 2. Stokes V LSD profiles of DT Vir=GJ 494A (thick black line) along with the maximum entropy fit (thin red line) to the data, 
for both the 2007 (left) and 2008 (right) runs. The rotational cycle of each observation, along with 3cr error bars, are shown next to each 
profile. 
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Figure 3. Magnetic topologies of DT Vir=GJ 494A in 2007 (top) and 2008 (bottom), reconstructed from a fit to the complete series of 
LSD Stokes V profiles. The three components of the field in spherical coordinates are displayed (from left to right), with magnetic fluxes 
labelled in G. The star is shown in flattened polar projection down to latitudes of —30°, with the equator depicted as a bold circle and 
parallels as dashed circles. Radial ticks around each plot indicate the phases of observation. 
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Table 7. Same as Table[2]for GJ 49. Rotational cycles are now computed using ephemeris HJD= 2,454, 300.0 + 18. GE". 



Date 


HJD 


UT 


^exp 


S/N 


C^LSD 


Cycle 




Vr 






(2,454,000+) 


(h:m:s) 


(s) 




(10-4/e) 




(G) 


(kms-1) 


2007 Jul 


27 


OUO.vJJ-C/UO 




4x900 


330 


2.1 


0.463 


—4.6 dz 4.1 


—6.00 


Jul 


28 


OUc/.UOtJZ J. 


03:11:41 


4x900 


450 


1.5 


0.518 


_i 9 n -1- 9 s 


—6.02 


Jul 


29 


O -LU. UO^O -L 


UO. -LU. 


4x900 


410 


1.6 


0.572 


_9n Q _|_ q 1 


—6.02 


Jul 


30 


311.63656 


03:13:25 


4x700 


370 


1.8 


0.626 


-28.4 ±3.5 


-6.01 


Jul 


31 


312.63031 


03:04:19 


4x700 


380 


1.7 


0.679 


-29.2 ±3.3 


-6.06 


Aug 


01 


313.63144 


03:05:51 


4x700 


300 


2.2 


0.733 


-21.7 ±4.3 


-6.01 


Aug 


03 


315.63927 


03:16:55 


4x700 


320 


2.5 


0.841 


-5.9 ±5.6 


-6.02 


Aug 


04 


316.63804 


03:15:03 


4x700 


380 


1.7 


0.894 


-1.6 ±3.4 


-6.04 


Aug 


05 


317.60746 


02:30:56 


4x700 


350 


1.9 


0.947 


2.9 ±3.7 


-6.01 


Aug 


09 


321.61210 


02:37:15 


4x700 


340 


2.0 


1.162 


-10.2 ±3.7 


-6.02 


Aug 


10 


322.63155 


03:05:10 


4x700 


330 


2.1 


1.217 


-10.6 ±3.9 


-6.02 


Aug 


11 


323.63437 


03:09:08 


4x700 


280 


2.5 


1.271 


-16.4 ±4.7 


-6.04 


Aug 


15 


327.62517 


02:55:34 


4x700 


350 


2.0 


1.485 


-10.3 ±3.8 


-6.03 


Aug 


18 


330.61820 


02:45:17 


4x700 


370 


1.8 


1.646 


-28.1 ±3.5 


-6.03 


Aug 


19 


331.57357 


01:40:56 


4x700 


340 


1.9 


1.698 


-20.2 ±3.6 


-6.06 


Aug 


31 


343.56027 


01:20:59 


4x700 


280 


2.4 


2.342 


-5.9 ±4.7 


-6.00 



each epoch respectively. The reconstructed magnetic maps 
are shown in Fig. [3] 

As obvious from Fig. [S] the magnetic topologies we de- 
rive contain a significant amount of toroidal field (62% and 
47% of the reconstructed magnetic energy in 2007 and 2008 
respectively); at both epochs, the reconstructed toroidal 
field shows up as a ring of counterclockwise field encircling 
the whole star. The poloidal field is more complex, especially 
in 2008 when more than 50% of the reconstructed poloidal 
field energy concentrates in orders with i > 3 (with dipole 
modes containing only 10%). The intrinsic evolution of the 
large-scale field topology between 2007 and 2008 is straight- 
forwardly visible in the image, especially on the poloidal 
field component whose spatial structure was much simpler 
in 2007 (64% of the poloidal field energy in dipole modes). 
At both epochs, the reconstructed poloidal field is mostly 
non axisymmetric (less than 20% of the energy concentrat- 
ing in m < £/2 modes). This information is summarised in 
Table [1 



5 DS LEO = GJ 410 = HD 95650 = HIP 53985 

DS Leo is a single MO dwarf with sharp spectral lines, 
located at an Hipparcos distance of 11.66 ±0.18 pc 
from the Sun Its RV is equal to -13.90 ± 0.10 kms"^ 
idever et al.l 



2002). Photometric variability was studied 



bv lFekel Henrv (2000 ) who detected cyclic variability at 
periods of 13.99 d and 15.71 d on different observing sea- 
sons and interpreted it as caused by stellar rotation (about 
14 d) coupled to surface differential rotation (modulating 
the observed photometric period as spots migrate to dif- 
ferent latitudes) . U sing the mass- luminosity relations of 
Delfosse et al. we estimate that the mass of DS Leo 

is 0.58 ± 0.02 Mo, i.e. very similar to that of DT Vir. The 
large amplitude rotational modulation that we observe for 
Stokes V profiles suggest that the inclination angle is not 
small; we therefore set i — 60° in the following imaging pro- 
cess. 

Stokes V data were collected at 2 epochs (same runs 



as for DT Vir, see Table [3]), providing again partial cov- 
erage of the rotation cycle at the first epoch but a dense 
and redundant coverage at the second epoch. Stokes V sig- 
natures are detected in almost all spectra, with longitudinal 
fields never exceeding strengths of 35 G. The RV we measure 
( — 14.05 ± 0.10 kms~^ at both epochs, with an internal dis- 
persion of^^03_^ms~^)isjii reasonbly good agreement with 
that of iNidever et al.l (|2002) . The rotational broadening of 
DS Leo is small, only slightly larger than that of GJ 205 
(having vsini = 1.0 — 1.5 kms~^, iReinersI l2007f ): using 
vsini — 2 ± 1 kms~^ provides a good fit to the Stokes / 
profiles and is compatible with the radius expected from 
theoret ical models (0.52 R©, see Table [1]), the rotation pe- 
riod of iFekel Henrv (2000) and the inclination angle we 
assumed (z = 60°). As for DT Vir, we do not need to use / 
for modelling the profiles of DS Leo. 

The solid-body-rotation period providing the best fit 
to the data is close to 14 d at both epochs; we there- 
fore used it to phase all spectra. As for DT Vir, the Bi 
curve in the 2008 data is showing apparently discrepant 
points for spectra collected at nearby phases but different 
cycles (e.g., Jan 03 and Feb 14, at rotation cycles 26.336 
and 29.328) as a likely result of the presence of surface 
differential rotation. This is confirmed by the fact that 
the full Stokes V 2008 data set cannot be fitted down to 
= 1 when assuming solid body rotation. Proceeding as 
for DT Vir, we obtain that Qeq = 0.465 ± 0.004 radd"^ 
and an = 0.076 ± 0.020 radd"^ at the surface of DS Leo, 
i.e., that the rotation periods at the equator and the pole are 
respectively equal t o 13.5 d and 16.1 d (b racketing both pho- 
tometric periods of lFekel k, Henrvll2000h . The photospheric 
shear of DS Leo is thus very similar to that of DT Vir, with 
the equator lapping the pole by one cycle every 83^20 ^• 

The optimal maximum entropy fit to the Stokes V data 
that we obtain (including the effect of differential rotation) 
is shown in Fig. |4]and corresponds to = 0-9, i.e., to a x? 
improvement over a non-magnetic model of xl5 and x8 for 
each epoch respectively. The reconstructed magnetic maps 
are shown in Fig. [5] 

The magnetic topologies we derive are predominantly 
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Figure 4. Same as Fig.[2]for DS Leo=GJ 410. 
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Figure 5. Same as Fig.[3]for DS Leo=GJ 410. 



toroidal (more than 80% of the reconstructed magnetic en- 
ergy at both epochs), with the reconstructed toroidal field 
showing up as a ring of clockwise field encircling the star. 
The poloidal field is much simpler than that of DT Vir 
(partly because of the lower spatial resolution resulting from 



the smaller v sin i) and consists mostly in a dipole (contain- 
ing more than 50% of the poloidal field energy) evolving 
from a mainly axisymmetric to a mainly non-axisymmetric 
configuration between 2007 and 2008. Event hough the frac- 
tional energy stored in the various field components re- 
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mained grossly stable (see Table [8]), the magnetic topol- 
ogy underwent significant temporal evolution between both 
epochs, e.g., with the toroidal field ring showing 2 extrema 
across the star in 2007 while it only shows one in 2008. 



6 CE BOO = GJ 569A = HIP 72944 

CE Boo is a young and active M2.5 dwarf with sharp spec- 
tral lines, located at an Hipparcos distance of 9.81 ±0.16 pc. 
It is the brightest member of a multiple (possibly quadru- 
ple) system with an estimated age of only 100 Myr 
(jSimon et all l20Q6l ). The companions are located about 
^" (i.e., 50 AU) away and consist of at least a brown 
dwarf binary and possi bly even a triple. Us ing the mass- 
luminosity relations of iDelfosse et aD (|2000h . we estimate 
a mass of 0. 48 =b 0.02 Mc:) The current RV is equal to 
-7.21 ± 0.10 (|Nidever et al.ll2002l ). From photometric vari- 



ability, "Kiraga & Stepien (2007') find that the rotation pe- 
riod is 13.7 d; this is surprisingly long for a star as young as 
CE Boo, even longer than the average period for M dwarfs 
of the young galactic di sc (whose age is typically a few Gyr, 
iKiraga &: SteDienll2007t ). The activity of CE Boo is slightly 
larger than that of DS Leo (but smaller than that of DT Vir, 
see Table [1]) , in agreement with what is expected for a star 
with a similar period and a later spectral type. The sharp 
lines of CE Boo also argue in favour of the slow rotation. 

Stokes V data were collected in early 2008 only, with 
a rather dense coverage of the rotation cycle (see Ta- 
ble 2]); Stokes V signatures are detected at all times, with 
longitudinal fields ranging from -50 G and -120 G and 
evolving smoothly with rotation phase. The RV we mea- 
sure (—7.34 ± 0.10 kms~^, with an internal d ispersion of 
0.03 kms"^) agrees with that of iNidever et al.1 (|2002l ). The 
rotational modulation of the Stokes V profiles are reminis- 
cent of those of AD Leo (M08), suggesting that the star 
is not seen equator on; the relative fluctuations of the lon- 
gitudinal fields are however about twice larger than those 
of AD Leo, indicating that i is not as low as 20° (as for 
AD Leo). We chose i — 45° as a an intermediate value. 
The rotational broadening in the spectral lines of CE Boo 
is small, comparable to that of DS Leo; using the radius 
expected from theore tical models (0.43 Rq , see Table [1]), 
the rotation period of iKiraga k, StepienI (2007 ) and the in- 
clination angle we assumed {i — 45°), we find (and used) 
vsini — 1 kms~^. Conversely to the 2 previous stars, we 
have to assume / = 0.05 (smaller than the usual value for 
mid-M dwarfs, see M08) to obtain a = 1 to the Stokes V 
profiles. 

The solid-body-rotation period providing the best fit to 
the Stokes V data is equal to 14.7 d, which we used to phase 
all our spectra. This is slightly longer than the period found 
bv lKiraga &: StepienI (|2007), suggesting that CE Boo is also 
subject to differential rotation (as DT Vir and DS Leo) with 
at least dO. — 0.03 radd~^ (assuming our period yields the 
rotation rate at the pole and the period of Kiraga & Ste pienI 
I2OO7I traces the rotation rate at the equator) . Proceeding as 
above, we obtain no clear minimum in the f^eq-<^^ domain, 
indicating that our data are not suitable for measuring dif- 
ferential rotation; this is not too surprising given the fairly 
simple rotational modulation of the Stokes V profiles and 
its moderate amplitude. We therefore assumed that CE Boo 
rotates as a solid-body in the following; very similar results 
are obtained if assuming that CE Boo is hosting differential 
rotation similar to that of DT Vir and DS Leo. 
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The optimal maximum entropy fit to the Stokes V data 
that we obtain is shown in Fig. [6] and corresponds to = 1, 
i.e., to a xi improvement over a non-magnetic model of x35. 
The reconstructed magnetic map is shown in Fig. [71 The 
magnetic topology we derive is almost completely poloidal, 



with less that 10% of the reconstructed energy concentrating 
into the toroidal component (see Table [8]); the poloidal field 
is quite simple (93% of the energy in the i < 2 modes) and 
mostly axisymmetric (96% of the energy in m < i/2 modes). 
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7 OT SER = GJ 9520 = HIP 75187 

OT Ser is an active Ml. 5 dwarf with spectral lines showing 
significant rotational broadening. Located at an Hipparcos 
distance of 11-39 di 0. 25 pc, it has no identified companion 
(| Paemgen et al ] l2007l . Forveille privat e communicati o n) . Us - 
ing the mass- luminosity relations of IPelfosse et al.l (|2000l ). 
we estimate a mass of 0.55 di 0.02 M©. Two discrepant ro- 
tation periods, both estimated from pho tometric variability, 
are reported i n the literature; while [ Norton et al.l (|2007l ) find 
Prot — 3.38 d. lKiraga Step ie 3 (|2007|) obtain Prot = 0.37 d. 
Given the estimated radius of OT Ser (about 0.5 R©, see Ta- 
ble [TJ, the second period would imply an equatorial velocity 
of almost 70 kms~^, much larger tha n the observed width 
of spectral lines; the rotation period of iNorton et al.l (|2007l ) 
is thus much more likely to be the correct one. 

Stokes V data were collected in 2007 and 2008, with a 
rather dense coverage of the rotation cycle in 2008 (see Ta- 
ble[5}; Stokes V signatures are detected in almost all spectra. 
The longitudinal field variations with phase are very differ- 
ent at both epochs (with Bi varying from 55 to 80 G in 2007 
and from -10 to 110 G in 2008) demonstrating that the mag- 
netic topology changed significantly on a timescale of only 
0.5 yr. The RV we measure (6.77 ± 0.10 kms"^ and 6.83 ± 
0.10 kms~^, with internal dispersions of 0.03 kms~^and 
0.04 kms~^) are slightly different at both epochs, possibly 
reflecting the change in the magnetic topology. Modelling 
Stokes / line profiles yields vsini = 6 it 1 kms~^. From the 
expected radius and rotation period, we infer that the star is 
seen at an intermediate inclination angle; we use i = 45° in 
the following. As for CE Boo, we have to adjust / to obtain 
a = 1 fit to the Stokes V profiles; we find that / equals 
0.05 and 0.10 in 2007 and 2008 respectively. 

The solid-body-rotation period providing the best fit to 
the data is equal to 3.40 d at both epochs; we used it to phase 
all our data (see Table [5]) . This period is cl ose to but slightly 
different than that of i Norton et al ] ([2007'), suggesting that 
OT Ser is also differentially rotating. Fitting our 2008 Stokes 
V data further confirms that OT Ser is not rotating as a solid 
body; with the same procedure as above, we obtain that 
^leq = 1-88±0.01 radd-^ and dO. = 0.12±0.02 radd-\ i.e., 
that the rotation periods at the equator and the pole are 
respectively equal to 3.34 d and 3.57 d. The photospheric 
shear of OT Ser is thus apparently even stronger than that 
of PT Vir and PS Leo, with the equator lapping the pole 
by one cycle every 521^^ d. 

The optimal maximum entropy fit to the Stokes V data 
that we obtain (including the effect of differential rotation) 
is shown in Fig. [8] and corresponds to = 1? i-^-, to a xi 
improvement over a non-magnetic model of x23 and xl9 
for 2007 and 2008 respectively. The reconstructed magnetic 
maps are shown in Fig. [9] Although both maps show a simi- 
lar large-scale topology (e.g., same latitudinal dependence of 
field polarities for all components), differences are neverthe- 
less obvious; for instance, the ring of positive radial field en- 
circling the star at mid latitudes shows a prominent blob at 
phase 0.0 in 2008 (causing the large- amplitude longitudinal- 
field modulation observed at this epoch). 

The magnetic topologies we derive are dominantly 
poloidal, with about 20-30% of the reconstructed energy 
concentrating into the toroidal component; the poloidal field 



is mostly axisymmetric and includes a significant dipole 
component at both epochs (see Table [8} . 



8 GJ 182 = HIP 23200 

GJ 182 is a very young single MO. 5 dwarf of the IC 2391 su- 
percluster located at an Hipparcos distance of 26.7 ± 1.7 pc. 
The star is surrounded by a mas sive debris disk indicat- 
ing on-going planetary formation (|Liu et al.ll2004l ). further 
demonstrating that it is indeed very young. Its position in 
the HR diagram (about 0.5 mag above the main sequence) is 
in agreement with the age of its yo ung moving group (about 
35 Myr, e .g.. iMontes et al. 200lh . Using the evolutionary 
models of iBaraffe et al.l ( 19981 ) and matching them to an 
absolute V magnitude and a logarithmic luminosity (rela- 
tive to the Sun) of 7.94 and —0.83 respectively, we find that 
GJ 182 has a mass of 0.75 M©, a radius of 0.82 i?^, a tem- 
perature of 3950 K and an age of 25 Myr; this is what we 
assume in the following. The high lithium content of GJ 182 
suggests that t he star is even yo unger, possibly as young 
as 10-15 Myr (|Favata et al.l Il998t ) as evolutionary models 
predict t hat lithium should be already strongly depleted at 
20 Myr (|Favata et al.lll998l ). Effects of rotation and mag- 
netic fields on the stellar structure and on the evolution 
(e.g., IChabrier et al. I I2OO7I . not taken into account in exist- 
ing studies) are however likely to affect model predictions 
significantly. 

Stokes V data were collected in 2007, covering only 
about half the rotation cycle of GJ 182 (see Table [6]); Stokes 
V signatures are detected in all spectra and longitudinal field 
vary from -90 G to 20 G with rotation phase. The RV we 
measure (19.35 db 0.10 kms~^, with an internal dispersion of 
0.18 kms~^) varies with rotational phase and correlate well 
with longitudinal field values; although the statistics is mod- 
erately significant (only 7 data points available) , it suggests 
that the RV fluctuations we detect (full amplitude of about 
0.4 kms"^) are related to surface magnetic activity. Spec- 
tral lines are significantly broadened by rotation; modelling 
Stokes / LSP profiles yield vsini = 10 ± 1 kms~^. Pho- 
tomet ric modulation indicates a rotation period of about 
4.4 d ('Kiraga k Stepienl [2OO7I ) . suggesting that the star is 
viewed equator-on rather than pole-on; we therefore set the 
inclination angle i — 60°. 

The solid-body-rotation period providing the best fit 
to the data is equal to 4.35 d, slightly smaller than the 
photometric period of 4.41 d measured bv lKiraga k Stepienl 
(2007); we used our estimate to phase all spectra, and take 
this as a likely indication that GJ 182 is a differential ro- 
tator. Pespite the small number of spectra and the limited 
phase coverage, fitting our Stokes V data down to x^ = 1 



suggests that GJ 182 is indeed not rotating as a solid body. 
Using the procedure described in Sec. [3l we obtain that 
Qeq = 1.46 ± 0.01 radd-^ and = 0.06 di 0.03 radd^^ 
the corresponding rotation periods at the equator and the 
pole are respectively equal to 4.30 d an d 4.49 d (bracket- 
ing both our rotation period and that of 'Kiraga k Ste pienl 
[2007). SoUd-body rotation is excluded at the 2a level; more 
data are needed to confirm this with better precision. 

The optimal maximum entropy fit to the Stokes V data 
that we obtain (including the effect of differential rotation) 
is shown in Fig. [TOland corresponds to xS = I5 i-e., to a xS 



Magnetic topologies of early M dwarfs 13 




n \ < \ 1 \ ^ 1 r 




-20 20 40 
Velocity (km/s) 
Figure 10. Same as Fig.[2]for GJ 182. 

improvement over a non-magnetic model of x22. The recon- 
structed magnetic map is shown in Fig. [11] The magnetic 
field is dominantly toroidal, and the poloidal component is 
mostly non-axisymmetric (see Table [5). 



9 GJ 49 = HIP 4872 

GJ 49 is a single Ml. 5 dwarf with sharp spectral lines and 
relatively low activity; it is the least active star of our 
sample (see Table [T]). The Hipparcos distance is equal to 
10.06 ± 0.14 pc. No rotation periods are fou nd in the liter- 
ature ; the RV is reported to be — 5.97± 0.10 (Nideve r et al. 
I2OO2I) . Using the mass- luminosity relations of Delfosse et al. 
(|200Q| ). we estimate a mass of 0.57 di 0.02 Mq. 

Stokes V data were collected in 2007, covering the whole 
rotation cycle of GJ 49 (see Table [7|); Stokes V signatures 
are detected in all spectra, with longitudinal fields rang- 



ing from -30 G to G across the cycle. The RV we mea- 
sure (—6.02 di 0.10 kms~^, with an internal dispersion of 
0.02 kms"^) is in good agreement with that of lNidever et all 
(2002 ). Modelhng Stokes / LSD profiles indicates that the 
rotational broadening is very small; we thus set v sin i = 
1 kms"\ 

We determine the rotation period by selecting the one 
with which the Stokes V profiles can be fitted at = 1 
with smallest magnetic energy in the reconstructed image; 
we find that Prot = 18.6it0.3 d; we also find that an inter- 
mediate inclination angle (i ^ 45°) minimises the amount of 
reconstructed magnetic information. We find that the data 
are compatible with solid-body rotation, but the accuracy 
to which we measure dQ (error bar ~ 0.05 radd~^) is not 
high enough to know whether GJ 49 also hosts differential 
rotation similar to that found on the other sample stars. 

The fit to the Stokes V data that we obtain is shown in 
Fig. [12] and corresponds to a xi improvement over a non- 
magnetic model of xll. The reconstructed magnetic map is 
shown in Fig. [13] The poloidal and toroidal field components 
roughly share the same amount of energy, with the poloidal 
field being mainly dipolar and axisymmetric (see Table [8}. 



10 SUMMARY AND DISCUSSION 

We report in this paper the results of our spectroscopic sur- 
vey of M dwarfs; following M08 (concentrating on mid-M 
dwarfs), we describe here the Zeeman signatures and the 
large-scale magnetic topologies we observed on 6 early-M 
dwarfs (from MO to M3). We also determined or confirmed 
the rotation period of all stars (ranging from 2.8 to 18.6 d), 
and detected significant surface differential rotation in 4 of 
them (with a strength comparable to that of the Sun). 

The magnetic fields we detect in early-M dwarfs are 
weak (typically a few tens of G), smaller in particular than 
those found in mid-M dwarfs (M08) by typically a factor of 5 
(see Fig. [13]) with a sharp transition occuring at 0.4 Mq (see 
Fig. 1151 left panel). The large-scale magnetic topologies we 
derive are also significantly different, involving a much larger 
fraction of toroidal fields and a lower axisymmetric degree 
of poloidal fields whenever > 0.5 Mq (5 stars in the 
present sample); below 0.5 Mq, the poloidal field is largely 
dominant and axisymmetric and its strength is increasing 
rapidly as mass decreases. We also observe that the typical 
lifetime of the large-scale magnetic topology is very different 
on both sides of the 0.4-0.5 Mq threshold, with lifetimes 
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smaller than a few months on the hot side and longer than 
1 yr on the cool side (M08). This threshold is very sharp 
and well defined, with little apparent dependence with the 
rotation period. 

At this stage, it is interesting to consider the effective 
Rossby number i?o, defined as Ro — Prot/'Tc where Tc is 
the convective turnover tim^B; in particular, Ro is a conve- 
nient parameter for comparing the strength of dynamo ac- 
tion (and, e.g., logi^x that indirectly reflects dynamo action 
through coronal heating) in stars with different masses. Fig- 
ure [15] (right panel) illustrates how log i?x smoothly varies 



^ The Tc values that we use here are those of i Kiraga StepienI 
determined empirically from relative X-ray luminosities 
of stars with different masses and rotation periods. At masses of 
~ 1 M0, they match the usual value of ~ 15 d; they steeply 
increase with decreasing mass below masses of ~ 0.6 M©. 



with Ro for both early- and mid-M dwarfs studied here and 
in M08; we find that log i?x increases steeply with decreas- 
ing Ro until Ro :^ 0.1 where logi?x saturates at a level of 
about -3.0, in go od agreement with previous studies (e.g., 
I James et al.]l2Q00l ). Only GJ 182 (at Ro = 0.17) lies slightly 
above the overall trend, as a likely consequence of its extreme 
youth and the related differences in its internal structure. 

We note that the abrupt step in the large-scale magnetic 
energy between early- and mid-M dwarfs (see Fig. I15[ left 
panel) apparently correlates better with stellar mass than 
with Ro] at Ro — 0.05, AD Leo exhibits a magnetic flux 
compatible with that of all other early M dwarfs, but sig- 
nificantly weaker than that of EV Lac (at Ro — 0.07) and 
YZ CMi (at Ro = 0.04). More data are needed to confirm 
this point, especially for high-mass rapid rotators (i.e., hav- 
ing Ro ^ 0.01 and > 0.5 M©) and low- mass slow ro- 
tators (with i?o > 0.1 and < 0.4 M©). We also note 
that this abrupt step does not show up in the log Rx vs Ro 
plot of Fig. [15] (right panel); this is presumably because X- 
rays are sensitive to overall magnetic energies while we are 
only sensitive to the largest scales. Our result therefore sug- 
gests that, at some specific stellar mass (^ 0.4 M©, rather 
than at some specific Ro), dynamo processes become sud- 
denly much more efficient at triggering large-scale magnetic 
fields; we also observe that, at more or less the same mass 
(~ 0.5 M0), large-scale topologies of M dwarfs become dom- 
inant ly poloidal and axisymmetric. Note however that, even 
in the case of mid-M dwarfs, the large-scale fields we derive 
are significantly smaller than the corresponding equiparti- 
tion field (a few kG, M08). 

Significant surface toroidal fields are detected even in 
DS Leo and GJ 49, i.e., the two slowest rotators with masses 
larger than 0.5 M©; it suggests that the transition between 
mainly poloidal and mainly toroidal fields in > 0.5 M© 
stars occurs at Ro :^ 0.5 — 1.0, with the Sun located on the 
other side of this boundary (at Ro :^ 1.5 — 2.0). Note that 
this boundary coincides with the sharp onset of photometric 
variability in convective stars (occuring below Ro 0.7, 
lHallll99lh . With a poloidal field concentrating 70-80% of 
the reconstructed magnetic energy, OT Ser is off this trend; 
we suspect that this is due to its proximity with the 0.5 M© 
sharp threshold below which magnetic topologies become 
dominantly poloidal. 

Early- M dwarfs are found to show significant differential 
rotation; the values we obtain for the surface angular rota- 
tion shear dQ ranges from 0.06 to 0.12 radd~^, i.e., from 
once to twice the strength of the surface latitudinal shear of 
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Figure 14. Basic properties of the large-scale magnetic properties of early- and mid-M dwarfs as a function of stellar mass and rotation 
period. Symbols size indicates magnetic energies, symbol colour illustrates the field configuration (blue and red for purely toroidal and 
purely poloidal fields respectively) while symbol shape depicts the degree of axisymmetry of the poloidal field component (decagon and 
stars for a purely axisymmetric and purely non-axisymmetric poloidal fields respectively). Results for early-M stars are from this paper 
and results for mid-M stars are from M08. The full, dashed and dash-dot lines respectively trace the location of the Ro = 1, 0.1 and 0.01 
contours, while the dotted line shows the theoretical full-convection threshold (M^ ~ 0.35 M©). 
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Figure 15. Reconstructed magnetic energy (left) and relative X-ray luminosity (with respect to the bolometric luminosity, right) as a 
function of Ro for stars studied in this paper and in M08. Stars with masses larger and smaller than 0.4 are shown as green squares 
and red circles respectively. In the left panel, measurements at different epochs (whenever available) are shown for each star to illustrate 
the typical scatter expected from temporal variability. 



the Sun. Our detection is further confirmed by the small 
(but significant) differences between the rotation periods 
we measure and the values reported in the literature (de- 
rived from photometric fluctuations) and by the short life- 
times of the large-scale field topologies (quickly distorted be- 
yond recognition by differential rotation). Previous Doppler 
imaging studies of early-M dwarfs with very fast rotation 



Ro :^ 0.01) report that differential rotation is very small 
Barnes et al.ll2005l l: our study suggests that the situation 
may significantly differ in moderate rotators like those we 
considered. Our result is also different from what is observed 
in mid-M dwarfs where differential rotation is very small 
(a few mradd"^ at most, i.e., more than 10 times smaller 
than that of early-M dwarfs, M08) and large-scale magnetic 
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Table 8. Properties of the large-scale field topologies derived in the present study. For each star, observations at different epochs are 
listed separately. The table lists sequentially the name of the star, the mass, the rotation period, the effective Rossby number and the 
logarithmic relative X-ray luminosity (taken from Table the surface angular rotation shear between the equator and pole dO. (whenever 
detected), the derived filling factor / (whenever applicable), the reconstructed magnetic energy and flux (i.e., < > and < B >), 
and the fractional energies in the poloidal fleld, the poloidal dipole (£ = 1) modes, the poloidal quadrupole (£ = 2) modes, the poloidal 
octupole (£ = 3) modes, and the poloidal axisymmetric (m < ^/2) modes. All fleld components with £ < 3 are required to flt the data at 
noise level. 
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GJ 182 2007 


0.75 


4.35 


0.174 


-3.1 


0.06 ±0.03 
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DS Leo 2007 
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0.20 


0.52 
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GJ 49 2007 
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0.09 
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0.48 
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OT Ser 2007 
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0.097 


-3.4 




0.05 


2.28 
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0.86 


2008 
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0.10 
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0.33 


0.17 
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CE Boo 2008 


0.48 


14.7 


0.350 


-3.7 




0.05 


1.48 
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0.95 


0.87 


0.06 
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0.96 



topologies long-lived (M08). It is not clear yet what this dif- 
ference is due to; while small Ro may contribute at freezing 
differential rotation, this is likely not the only relevant pa- 
rameter for this problem (e.g., with DT Vir and EV Lac 
showing respectively significant and no differential rotation 
despite their similar Ro). 

The sharp transition that we report between the mag- 
netic (and differential rotation) properties of early- and mid- 
M dwarfs is surprising at first glance; naively, one would ex- 
pect the properties of large-scale magnetic fields to change 
smoothly with stellar mass as the radiative core ^ ets progres- 
sively smaller. From the evolutionary models of ISiess et alJ 
([2OOQ), we however note that the outer radius of the radia- 
tive core of early- M dwarfs is changing very quickly with 
stellar mass, from about 0.5 Ri, for a 0.5 Mq star to a neg- 
ligible fraction for a 0.4 Mq star. We speculate that this 
sharp transition is the main reason for the abrupt magnetic 
threshold that we report here. The rapid increase in empir- 
ical conv ective turnover times oc curing at about the same 
location teiraga SteDienl l2007l l also likely contributes at 
making the transition between both dynamo regimes very 
sharp. 

The most recent dynamo simulations of fully convec- 
tive M dwarfs (Browning 2008) (carried out for Ro ^ 0.01) 
are successful at reproducing the frozen differential rotation 
that we observe (MOB); they however predict the presence 
of strong toroidal fields that we do not see in mid-M dwarfs 
with similar Ro. We speculate that the abrupt change in 
the large-scale magnetic topology of M dwarfs that we re- 
port here to occur at spectral type M3 may also be (at least 
partly) responsible for the reduced magnetic braking ob - 
served for stars later than M3 (e.g., IPelfosse et al.lll998h : 
MHD simulations of magnetic winds are necessary to es- 
timate quantitatively whether the observed change in the 
large-scale magnetic topology can indeed explain the longer 
spin-down timescales. 

In the two stars having v sin i measured with sufficient 
precision (i.e., vsini > 10, GJ 182 and DT Vir), we find 
that R^sini (equal to 0.86 ±0.09 R© and 0.62 ±0.06 R© re- 
spectively) is alrea dy larger than the p redicted radius from 
theoretical models (jBaraffe et al.|[l998h : while this could re- 



sult from overestimating the true age (and hence underesti- 
mating the true radius) of GJ 182, this explanation does not 
apply for DT Vir, for which we conclude that the observed 
radius is truly larger (by at least 10% and potentially as 
much as 30%) that what theoretical models predict. A simi- 
lar conclusion is re ached for V374 Peg (|Donati et al.ll2006al : 
iMorin et al1l2008al ): furthermore, MOB obtains that R^, sin i 
is equal to the predicted theoretical radius (within the error 
bars) for 4 other active stars, suggesting again that R j, is 
larger than expected. Following Chabrie r et al.l (|2007f ), we 
propose that this effect is a direct consequence of magnetic 
fields getting strong enough (and hence saturating the dy- 
namo, see Fig. [15]) to affect the energy transport throughout 
the convective zone and hence the radius. Our results there- 
fore independently confi rm the report that c ool low-mass ac- 
tive stars, either s ingle dMorales et al.ll2008l l or within close 
ecUpsing binaries (|Ribasi]2"o06l V usuallv have oversized radii 
with respect to inactive stars of similar spectral types. 

We also detect significant RV fluctuations (with a full 
amplitude of up to 0.40 kms~^) in the 3 very active stars of 
our sample (with rotation periods smaller than 5 d) . For the 
most active ones (DT Vir and GJ 182, showing the largest 
RV modulation) , the RV variations correlate reasonably well 
(though not perfectly) with longitudinal fields, suggesting 
that the origin of the variations is indeed the magnetic field 
(and the underlying activity). It also suggests that spec- 
tropolarimetric observations should be carried out simulta- 
neously with RV measurements of active stars to enable fil- 
tering out efficiently the activity jitter from the RV signal; 
this technique may prove especially useful when looking at 
Earth-like habitable planets orbiting around M dwarfs in the 
future, e.g., with a nIR spectropolarimeter such as SPIRou 
(a nIR counterpart of ESPaDOnS, proposed for CFHT). 

Our spectropolarimetric survey is an on-going study; we 
are now concentrating on late-M dwarfs (M5-M8) to derive 
similar observational constraints about the large-scale mag- 
netic topologies of stars in the yet unexplored 0.08-0.20 M© 
region of Fig. [T3] to investigate how dynamo processes op- 
erate down to the brown dwarf threshold, i.e., when stellar 
atmospheres get so cool that they start to decouple from 
their magnetic fields. 
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